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Executive  Summary 

In  two  previous  reports(i'2),  a  proposal  was  tabled  for  designing 
and  certifying  aerospace  structures  of  composite  laminate,  and  for 
repairing  cracked  aluminum  structures  using  composite  patches.  In 
particular,  it  was  based  on  sound  fundamental  understanding  of  the 
mechanisms  by  which  damage  accumulates  in  carbon  fiber  composites 
and  structure  across  several  orders  of  magnitude  of  size:  fiber  breakage 
and  matrix  cracking  that  lead  progressively  to  ply  failure;  delamination 
and  laminate  cracking  that  lead  to  component  fracture.  Application  of 
this  approach  to  predictive  design  requires  the  necessary  "know-how" 
and  procedure  for  dealing  with  die  identification  of  failure  mechanisms 
from  the  micron  to  the  metre  level  of  size. 

One  method  which  is  empirical  in  nature  is  to  assess  in  sequence  the 
effect  of  cyclic  stress  on  the  residual  strength  and  lifetime  of  a  test 
coupon,  substructural  element,  component  and  finally  large  scale 
structure  by  a  progression  of  experiments.  This  methodology  is  based  on 
transferring  empirical  information  and  relating  experimental  data  from 
one  point  on  the  size  scale  to  another.  Empirical  methods,  however,  do 
have  serious  limitations  and  shortcomings. 

The  proposed  way  toward  is  to  begin  the  design  process  at  the 
constituent  level  with  physical  model-informed  selection  of  fiber  and  it's 
surface  treatment,  matrix  and  it's  processing  condition,  and  choice  of 
preferred  architecture  (geometry)  of  laminate  construction.  As  the  shift 
of  information  bothways  along  the  size  scale  proceeds,  the  design 
process  at  each  level  of  size  evolves  and  includes  the  identification  of 
the  principal  failure  mechanisms,  delamination,  fiber  fracture,  and  so 
on.  A  failure  mechanism  missed  or  not  understood  can  have  disastrous 
consequences  on  the  life  of  the  structure. 

Making  links  or  connections  between  our  understanding  of  material 
and  component  behaviour  relies  critically  on  this  knowledge  of 
microscopic  and  macroscopic  structural  changes  over  the  entire  span  of 
size  and  lifetime  of  the  part.  Hence,  a  technique  for  the  evaluation  of  the 
conditions  for  growth  (or  no  growth)  of  matrix  cracks,  delaminations 
and  fiber  fractures,  for  instance,  is  an  essential  tool  in  the  determination 
of  damage  tolerance  of  composite  material  and  structure.  Numerical 
methods  are  generally  available  through  the  use  of  finite  element 
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analysis  and  the  crack  closure  integral  technique  based  on  fracture 
mechanics  principles. 

Whilst  a  fracture  mechanics  approach  is  not 

mature  at  present  to  warrant  a  recommendation  for  genera 
application  to  the  certification  of  developmental  composite  hardware, 
nSrss  tWs  report  considers  the  merits  of  fracture  mech^cs 
Principles  with  respect  to  the  important  failure  mechanism  of  fibe 

Lcture,  including  fiber  breaks  precipitated  ^  be“ttIat^’usiM  t^ 
cracks  Individual  cracks  of  these  kinds  can  be  treated  using  m 
principles  of  fracture  mechanics;  collectively,  their  effecte  on  the 
properties  of  the  material  and  strucmre  can  be  quantified  usmg 
ielafcnships  evolved  from  damage  mechanics  ideasd-^). 

Tn  constructing  a  physical  model,  use  is  made  of  the  tools  of  engineermg 
fn^Ss  science:  the  equations  and  principles  of  mechamc^ 
MdS^c^  kinetics,  and  so  on.  When  the  model  is  combmed  with 
*perim^tal  data  of  (say)  fatigue  strength, 

design  methodology  having  the  power  of  prediction  of  lifetime  wnic 
S  fiom  an  understanding  of  the  relevant  (rneanmg  dommant) 
mechanism(s).  At  best,  the  model  ^"capsulates  toe  ph^ic^  bas^s  o^h 
problem  in  mathematical  form;  it  summaries  '^^at  has  been  obse 
and  predicts  behaviour  under  conditions  that  have 

investigated. 

The  need  to  model  the  dominant  mechanism(s)  and  to  incorporate  the 
rnodd  in  toe  design  process  itself  need  only  be  pertment  m  the  fet 
McnaTtoe  "itocrLcopic"  size  level.  It  is  important  to  consider 
however,  that  the  local  tensUe  strength  of  the  0°  P^^ 
zone  which  ultimately  determines  lammate  performance  y 

affected  b^^  presence  of  matrix-dominated  ply  cracks  in  the  damage 

(2)  the  presence  of  a  hostile  environment  which  may  lead  to  stress- 

“"IsTcTaS^s'^rtemperature  that  set  up  thermal  and  residu^ 
stressL  For  example,  there  is  convincing  evidence  widdy 
that  the  strength  of  glass  fiber  depends  on  the  concentration  of 
M  thSoXrrLd  the  degree  of  acidity  or  alkalinity.  Also,  that  toe 

shaje  of  damage  zones  in  carbon  fiber-polymer 

to  temperature  which  affects  matrix  ductility,  toughness,  and  so  forth. 
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In  reports  1  and  2, 1  considered  matrix-dominated  failure  processes  that 
can  occur  predominantly  in  carbon  fiber-epoxy  laminates  (CFRP)^ 
essentially  splitting  and  delamination.  In  this  Third  Report,  I  consider 
an  additionl  mechanism  of  failure  that  is  dominant  in  glass  fiber-epoxy 
laminates  (GFRP)  that  of  fiber  fracture.  In  particular,  this  Third  Report 
proposes  a  new  physical-based  micromechanical  theory  of  stress- 
corrosion  cracking  in  aligned  continuous  glass  fiber-epoxy  laminates. 

This  new  theory  is  based  on  the  premise  that  a  crack  initiates  at  a  pre¬ 
existing  surface  flaw  in  glass  fiber  and  extends  by  environmentally 
assisted  slow  crack  growth.  A  prediction  of  macroscopic  crack  growth  in 
the  laminate  by  simultaneous  matrix  failure  and  fiber  breakage  is 
derived  theoretically  as  a  function  of  the  crack  tip  stress  intensity  factor. 
Theoretical  prediction  which  compares  favourably  with  experimental 
measurement  indicates  that  crack  growth  rate  depends  on  the  size  of 
inherent  surface  flaw  of  glass  fiber,  the  fiber-matrix  interfacial  bond 
strength,  and  matrix  toughness. 

Finally,  for  completeness,  I  have  presented  the  essential  ingredients  of 
the  splitting  and  delamination  models(^'2)  and  equations  for  CFRP  as  an 
Appendix  to  this  Report. 
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1.  ABrief  Overview  of  Design  Philosophy 


Traditionally,  the  process  of  design  has  involved  specialised 
treatment  of  fatigue  and  corrosion  fatigue  phenomena.  For  metals  this 
has  meant  obeying  Miners  rule  and  making  good  use  of  Goodinan 
diagrams.  For  fatigue  crack  growth  and  stress-corrosion  crackmg,  t  e 
Paris  law  of  crack  propagation  has  been  apphed  with  marked  success. 
In  sharp  contrast,  the  treatment  of  sinular  phenomena  for  advanced 
composite  structures  is  more  complicated  becaiise  ^ 

formation  of  competing  and  interacting  cracks  m  the  ^atenaL  (Meta^ 
generally  fail  by  the  growth  of  a  single  crack).  Frequently,  this  has  led  to 
composites  industry  experiencing  hardware  development  problen^ 
due  to  design  oversights  that  lead  to  a  multitude  of  matrix  cracks,  sphts 
and  delaminations,  and  fiber  fractures. 


Furthermore,  fundamental  differences  in  the  elastic  properties  and 
fracture  toughness  between  fiber  and  matrix  constituents  preven 
meaningful  comparisons  between  material  systems  behaviour  and 
damage  tolerance  although  the  nature  of  the  fiber /matrix  bond  appears 
to  play  a  significant  role  in  the  cracking  and  fracture  processes. 
Consequently,  the  characteristics  of  damage  accumulation  from  ^e 
material  system  to  the  next  are  not  easily  predictable  due  to  this 
multiplicity  and  complexity  of  interacting  cracks  and  competmg  failure 

modes. 


An  alternative  innovative  way  foward  is  to  begin  the  design  process  at 
the  constituent  level  with  physical  model-informed  selection  of  fiber  and 
it's  surface  treatment,  matrix  and  it's  processing  conditions,  and 
architecture  (geometry)  of  laminate  construction^'  ^  ^  j 

information  both  ways  along  a  size  scale  proceeds,  from  mdividual  fiber 
to  ply  and  from  laminate  to  structure,  the  design  process  at  each  level  of 
size  has  to  include  the  identification  of  the  principal  failure  mechamsms. 

Where  to  begin  the  modelling  process  on  the  size  scale  depends  np^^r  ^e 
usefulness  of  that  model  and  what  it  intends  to  accomplish.  Makmg 
links  or  connections  between  our  understanding  of  material  and 
component  behaviour  at  the  various  scale  of  sizes  relies  critical  y  on 
information  flow  in  both  directions  between  higher  and  lower  levels  of 
size,  and  a  working  knowledge  of  microscopic 
structural  changes  over  the  entire  span  of  size  and  lifetime  of  me 

component. 
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All  to  frequently,  the  effect  of  local  details  are  identified  very  late  in  the 
design  cycle  and  unfortunate  experiences  of  this  kind  contribute  in  large 
part  to  early  over  optimism  in  the  material's  initial  performance  figures 
and  design  rating.  Furthermore,  any  structural  analysis  performed  at 
this  stage  has  rarely  been  made  in  sufficient  detail  to  adequately 
evaluate  the  effect  of  potential  failure  modes  and  margins-of-safety. 
Thus,  those  failure  processes  influenced  by,  for  example,  matrix 
properties  which  have  a  strong  dependence  on  processing  conditions 
and  environment  are  susceptable  to  large  statistical  scatter  and  to  cyclic 
loading  effects. 

Potential  damage  initiation  sites  include  free  edges,  bolt  holes,  ply 
terminations,  and,  with  respect  to  repair  of  structures,  at  a  composite 
patch  interface  and  patch  itself.  Hence,  a  technique  for  the  evaluation  of 
the  conditions  for  growth  (or  no  growth)  of  cracks  is  an  essential  tool 
for  the  determination  of  the  damage  tolerance  of  composite  material 
and  structure.  Numerical  methods  are  generally  available  through  the 
use  of  finite  element  analysis  and  the  crack  closure  integral  technique 
based  on  fracture  mechanics  principles. 

The  application  of  fracture  mechanics  requires  a  significant  analysis  and 
test  effort  to  evaluate  "hot  spots"  within  the  structure  and  to  generate 
the  necessary  fracture  toughness  data.  Whilst  a  fracture  mechanics 
approach  is,  therefore,  not  considered  sufficiently  mature  at  present  to 
warrant  a  recommendation  for  general  application  to  the  certification 
of  developmental  composite  hardware,  nevertheless,  I  wiU  consider  the 
merits  of  fracture  mechanics  and  damage  mechanics  principles. 


2.  A  Brief  Consideration  of  Fracture  and  Damage  Mechanics 

In  most  advanced  composites,  the  earliest  signs  of  damage 
comprises  an  array  of  matrix  cracks  in  the  individual  transverse  plies, 
often  referred  to  as  first  ply  failure.  Following  transverse  ply  cracking, 
delamination  cracks  accompanied  by  splits  frequently  form,  including 
fiber  breaks,  precipitated  by  these  transverse  ply  cracks(i/  In  reality, 
transverse  ply  cracks  may  already  exist  in  the  virgin  material  caused  by 
a  combination  of  thermal  effects  during  processing  and  cooling  and  the 
mismatch  of  thermal  expansion  coefficients  between  adjacent  layers  of 
the  laminate.  Individual  cracks  of  these  kinds  can  be  treated  using  the 
principles  of  fracture  mechanics;  collectively,  their  effects  on  the 


-8- 


Towards  a  Predictive  Design  Methodology  of  Composite  Laminates  and  Patches 


properties  of  the  material  can  be  quantified  using  relationships  based  on 
damage  mechanics  ideas(i'  2). 

Damage  and  fracture  mechanics  principles  are  based  on  quantitative 
relationships  between  microscopic  and  macrosopic  parameters  that 
attempt  to  "bridge"  a  gap  that  has  opened  up  on  the  size  scak  between 
classical  materials  science  that  focuses,  for  example,  on  single  fiber 
reinforcement  models  at  the  microstructural  level,  and  structural 
engineering  and  hardware  development  at  the  one  metre  level  arid 
above.  It  is  precisely  at  this  gap  that  "design"  normally  changes  haiids 
from  being  a  materials  design  issue  to  being  a  structural  one  for  which 
neither  group  has  taken  responsibility. 

In  the  stictest  sense,  "damage  mechanics"  refers  to  the  physical 
mechanisms  of  multiple  cracking  by  which  damage  fori^  and 
propagates  in  the  material.  The  concept  of  continuum  damage 
mechanics  is  a  method  of  treating  a  "whole  population  of  microscopic 
defects".  In  expressing  that  damage  mechanics  and  fracture  mechaiucs 
are  complementary,  one  can  consider  that  the  final  stage  of  damage 
mechanics  is  the  initial  stage  of  fracture  mechanics  where  single  crack 
growth  of  a  major  defective  region  or  volume  of  material  takes  over. 

Turning  briefly  at  this  particular  point  in  the  report  to  a  related  and  also 
key  aspect  of  design  is  the  manner  in  which  a  material  under  stress 
reacts  to  the  presence  of  stress  concentrators  such  as  holes  or  notches. 
Whilst  a  number  of  macroscopic  empirical  models(i'  2)  offer  the  potential 
of  working  models  to  account  for  notch  size  effect  they  are  considerably 
restricted  in  their  application;  for  example,  they  do  not  provide  insight 
into  the  mechanisms  or  causes  of  damage  or  the  mechanics  of  ultimate 
failure.  Neither  do  they  predict  a  fracture  stress  for  different  intrinsic 
material  variables  or  extrinsic  (environmental  and  thermal)  conditions 
other  than  those  of  the  test;  or  provide  insight  into  the  connection 
between  data  at  the  laminate  and  microstructural  levels;  or  even  the 
reasoning  behind  the  selection  of  fiber  and  matrix  and  process  route  in 
the  first  place.  Such  a  design  approach  imposes  a  heavy  burden  of 
testing  and  data  collection  which  is  time-consuming  and  uneconomic  in 
this  current  cost-conscious  climate. 

A  logical  way  to  progress  might  be  to  express  the  material  ^d  critical 
stress  parameters  of  such  correlations,  if  possible,  as  functions  of  the 
physical  parameters  that  dictate  failure.  Some  of  these  baseline 
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parameters  can  be  described  using  the  alternative  approach  to  which  I 
have  drawn  attention  before(i'  2)^  that  based  on  the  physical  modelling 
of  actual  failure  processes  and  the  effects  these  mechanisms  have  on 
localised  stress  in  neighbouring  load  bearing  plies.  Interaction  processes 
between  mechanisms  become  extremely  important  since  they  determine 
the  strength  and  fatigue  lifetime  of  material  and  structure.  They  can 
best  be  understood  if  a  physical  model  of  each  dominant  mechanism  is 
constructed^!'  2), 


3.  The  Basis  of  a  Model 

A  model  is  a  gross  simplification  of  the  mechanism  but 
containing  the  essential  elements  of  which  there  are  many  types(!'  2): 

(1)  empirical  models  which  provide  a  mathematical  description  to 
a  set  of  experimental  data,  with  no  power  of  prediction  (an  equation  to 
fit  the  S-N  curve  of  fracture  stress  versus  load  cycles  to  failure  is  a  good 
example); 

(2)  physical  models  which  have  a  predictive  capacity  since  they 
are  based  on  well-known  physical  processes  (the  mechanisms  of 
delamination,  splitting,  and  fiber  fracture  within  a  laminate  are 
examples);  and 

(3)  static  and  dynamic  models  which  contain  material  properties 
that  remain  constant  with  time  (the  former)  or  show  time  or  rate 
dependence  (the  latter).  An  example  of  the  dynamic  model  is  the  rate- 
dependent  behaviour  of  Kevlar  fibers  in  epoxy  due  to  a  kinetic  process 
of  molecular  imkinking  in  the  fiber  imder  mono  tonic  or  cyclic  loading. 
Another  example  is  the  one  described  in  this  report,  namely,  stress- 
corrosion  (or  chemically  assisted)  cracking  of  glass  fiber-epoxy 
laminates. 

At  the  heart  of  the  model,  then,  lies  the  physical  mechanism(s)  which  can 
best  be  identified  by  direct  observation  using,  for  instance,  optical  or  in- 
situ  scanning  electron  microscopy.  It  is  dangerous  to  assume  a 
mechanism  without  evidence  that  it  actually  operates  and  is,  or  can  be, 
dominant  (meaning  that  it  has  more  influence  than  other  competing 
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In  constucting  a  model,  use  is  made  of  the  tools  of  engineering  and 
materials  science:  the  equations  and  principles  of  mechanics, 
thermodynamics,  kinetics,  and  so  on.  And,  of  course,  there  is  the  menu 
of  previously  modelled  problems  of  fracture  (or  damage)  mechamcs,  of 
thermal  stress  and  chemical  (envirorunental)  effects,  and  so  forth.  Many 
of  the  models  have  an  empirical  component  which  requires  a  more 
subtle  approach  or  "fine  tuning".  An  element  of  empiricism  introduces 
unknown  constants  into  the  equation.  Often,  they  simply  muhiply  to 
produce  only  one,  perhaps  two,  unknown  constants.  A  calibration 
procedure  is  then  required  to  set  values  on  them;  one  experimental 
measurement  for  each  constant.  Fine  tuning,  then,  simply  means 
comparing  the  model's  prediction  to  experimental  measurements  in  a 
significant  way,  covering  as  wide  a  range  of  conditions  as  possible,  fiber 
and  matrix  types,  fiber  lay-up,  stress-state,  number  of  fatigue  cycles, 
temperature,  environment,  and  so  on. 

The  design  of  critical  experiments  is  a  vital  element  to  this  process.  We 
can  never  get  away  from  doing  the  pertinent  experiment  which  is  often 
pointed  out  and  defined  for  us  by  fiie  model  itself.  When  the  model  is 
combined  with  experimental  data  of  (say)  fatigue  strength,  this 
approach  leads  to  a  design  methodology  having  the  power  of  prediction 
which  comes  from  an  understanding  of  the  relevant  (meaning 
dominant)  mechanims. 

At  best,  the  model  encapsulates  the  physical  basis  of  the  problem  in 
mathematical  form;  it  summaries  what  has  been  observed  and  predicts 
behaviour  under  conditions  that  have  not  been  investigated.  Caution  is 
imperative  when  extrapolating  the  model;  an  unknown  mechanism  not 
previously  encountered  can  have  disastrous  consequences.  Microscopic 
models  can  provide  insight  into  mechanisms  of  failure  and  in  parallel 
with  empirical  methods  they  can  be  used  to  develop  equations  having 
the  predictive  power  of  physical  modelling  with  the  precision  of  curve¬ 
fitting.  This  approach  is  given  respectability  by  having  been  assigned  the 
name  "model-informed  empiricism"  and  it  is  this  which  has  led  to  the 
development  of  the  new  branch  of  mechanics  called  "damage 
mechanics" . 

The  physical  modelling  approach,  then,  enables  us  to  gain 
understanding  of  the  origin  of  failure  processes;  to  capture  the 
material's  response  in  an  equation  or  code  of  design;  to  predict  material 
behavior  under  conditions  not  easily  simulated  in  the  laboratory  (for 
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example,  extremely  long  duration  fatigue  tests);  and  to  guide  the 
optimisation  of  micro(macro)structure  of  the  material  for  longevity  and 
durability.  The  aim,  then,  of  this  approach  is  to  establish  a  physical 
framework  within  which  empirical  descriptions  of  the  behaviour  of 
some  (preferably  all)  of  the  experimental  variables  in  composite  design 
are  embedded. 

A  word  of  warning:  such  procedures  of  developing  a  physical  model 
(with  empirical  ingredients,  as  appropriate),  may  not  work  equally  well 
for  all  systems.  We  are  at  a  starting  point.  For  comparison  with  the 
physical  models  of  matrix-dominated  failure  processes  in  carbon  fiber- 
epoxy  laminates(i'  2),  in  this  report  I  will  address  the  development  of  a 
new  micromechanical  theory  of  stress-corrosion  cracking  of  glass  fiber- 
epoxy  laminates  based  on  the  modelling  of  fiber  fracture. 


4.  A  Physical  Model  of  Fiber  Fracture  Based  on  Direct  Observation 

Glass  fiber-epoxy  laminates  can  be  weakened  by  the 
combined  influence  of  stress  and  corrosive  environment.  For  example, 
their  susceptibility  to  stress-corrosion  cracking  in  acidic  environments 
has  received  wide  attention(3-i7).  Hogg  and  HulK^),  Noble  et  al.(®)  and 
Price  and  Hull(9d2)  have  all  observed  stress-corrosion  cracking  of 
individual  glass  fibers  in  laminates  using  scanning  electron  microscopy. 
They  observed  that  a  micron-sized  crack  in  glass  fiber  initiates  at  a 
surface  flaw  and  propagates  stably  to  begin  with,  followed  by  unstable 
(fast)  propagation  across  the  remainder  of  the  fiber's  diameter.  The 
crack  then  extends  further  into  the  laminate  as  a  macroscopic  crack  that 
can  be  seen  using  an  optical  microscope.  In  fracture  mechanics  terms, 
the  relationship  between  macroscopic  crack  propagation  rate,  da/ dt, 
and  stress  intensity  factor,  Kl,  at  the  crack  tip  can  be  expressed  by: 

—  =  (1) 
dt  ' 

where  A  and  m  are  constants  which  depend  on  the  material  and 
environment. 

For  a  slowly  propagating  macroscopic  crack  there  is  little  or  no  fiber 
pull-out.  The  fracture  planes  of  the  fiber  and  matrix  are  coincidental. 
With  increasing  crack  propagation  rate,  there  is  an  increasing  amoimt 
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of  fiber-matrix  interfacial  debonding.  As  the  macroscopic  crack 
propagation  rate  increases  this  is  accompanied  by  an  increase  in  fiber- 
matrix  interfacial  debonding  indicated  by  longer  fiber  pull-out. 

The  formation  of  matrix  bridges  in  the  crack  wake  are  observed  when 
tougher,  more  ductile  matrices  are  used.  Such  mechanisms  that  act  in 
the  crack  wake  shield  the  tip  of  the  crack  from  localised  tensile  stress. 

In  the  next  section,  a  theoretical  relationship  between  the  macroscopic 
crack  propagation  rate  and  the  crack  tip  stress  intensity  factor  is 
deduced.  Comparison  is  made  between  theoretical  prediction  ^d 
experimental  data  which  includes  the  effects  of  fiber-matrix  debonding 
at  the  crack  tip  and  matrix  bridging  in  the  crack  wake. 


5.  A  New  Micromechanical  Theory  of  Crack  Propagation  in  GFRP 
(following  Sekine  and  Beaumont,  Cambridge  (1998),  unpublished) 

Following  the  work  of  Wiederhom  and  Bolz(i^),  the  crack 
propagation  rate,  da/dt,  in  bulk  glass  due  to  stress-corrosion  cracking 
is: 


da  (  E-ocKi 


(2) 


E  is  the  activation  energy  of  the  chemically-activated  process,  Ki  is  the 
stress  intensity  factor,  R  is  the  gas  constant,  T  is  absolute  temperature, 
and  V  and  Of  are  empirical  constants. 

On  the  surface  of  glass  fiber  are  large  numbers  of  flaws  of  variable 
size(io,  21).  It  is  reasonable  to  presume  that  a  crack  initiates  at  one  of 
these  pre-existing  surface  flaws  and  ^ows  perpendicular  to  the  fiber  s 
length  with  time  which  leads  to  a  time-dependence  of  fiber  strength 

(Fig- 1  )• 

Figure  2  shows  a  micrograph  of  a  typical  laminate  fracture  surface 
obtained  by  scanning  electron  microscopy.  The  fracture  surface  of  each 
glass  fiber  is  characterised  by  a  smooth  or  "mirror"  region  which 
indicates  slow  cracking  which  is  surrounded  by  a  "hackle"  region  which 
indicates  fast  fracture.  In  the  model,  the  shape  of  the  slowly  extending 
crack  front  is  represented  as  a  circular  arc  of  radius  r,  where  r  is  equal  to 
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the  fiber  radius.  The  average  crack  propagation  rate  due  to  some  stress- 
corrosion  process  can  thus  be  written  (from  eqn  (2))  as: 


r  1  Ydx\ 

v2r0jv  dt  j 


vexp 


E-aKi'" 
RT  y 


(3) 


where  X  is  the  imcracked  area  of  the  fiber,  0  is  half  the  angle  which  is 
made  by  two  fiber  radii  on  the  edges  of  the  stress-corrosion  crack,  and  t 
is  time  (Fig.  3).  In  eqn  (3),  the  stress  intensity  factor  Ki,  should  be 
interpreted  as  the  average  value  of  Ki  along  the  entire  length  of  the 
stress-corrosion  crack  tip.  Since  the  stress  intensity  factor  for  the  crack 
opening  mode  (designated  by  the  subscript  i)  is  constant,  more  or  less, 
along  the  entire  circular  crack  front,  Ki  can  be  represented  by  the  stress 
intensity  factor  at  the  maximum  depth  of  the  stress-corrosion  crack(22): 


Ki  =  aF{Q)^nKr  (4) 

where  (7  is  the  tensile  stress  acting  on  the  fiber  and  F(  0)  is  written  as 

F(0)=  V1-COS0  {l.l2-3.40(l-cos  0)+13.87(l-cos  0)2-14.37(l-cos  6)^}  (5) 

For  a  model  in  which  the  aligned  continuous  glass  fiber  are  distributed 
in  a  doubly  periodic  square  array  and  perpendicular  to  the  macroscopic 
stress-corrosion  crack,  the  relationship  between  the  tensile  stress  and 
the  stress  intensity  factor  for  the  crack  opening  mode  is  given  by(i2); 


<y=PK, 


(6) 


where: 


1  2vf 


(7) 


Here  Vf  is  the  volume  fraction  of  glass  fiber  and  V  is  that  fraction  of 
plies  of  the  laminate  in  which  glass  fiber  are  perpendicular  to  the 
macroscopic  crack. 
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Furthermore,  geometrical  consideration  of  the  area  of  the  stress- 
corrosion  crack  (shown  in  Fig.  3)  gives: 

^=— =-4r2sin2  6»—  (8) 

dt  dt  dt 

where  Y  is  the  area  of  the  stress-corrosion  crack  of  the  glass  fiber. 

Substituting  eqn  (8)  into  eqn  (3)  and  combining  eqns  (4)  and  (6),  we 
obtain: 


where: 


vk  6 


aj3F(^)V2;rr 


RT 


dO 


) 


k=QXp\ 


f 


— 1 
KT  ^ 


(9) 


(10) 


By  integrating  eqn  (9),  the  time  required  for  the  slow  crack  ^owth  stage 
of  failure  of  a  single  glass  fiber  by  stress-corrosion  cracking  is  given  by: 


f 


a^F{6)^27tr 

FT 


dO 


J 


(11) 


where  is  half  the  angle  made  by  two  fiber  radii  on  the  edges  of  the 
pre-existing  surface  flaw  and  6^  is  the  critical  value  at  unstable  (fast) 
fracture  of  the  glass  fiber. 

Hogg  and  Hull(3),  Aveston  and  Sillwood(7)  pointed  out  that  a  feature  of 
a  slowly  moving  macroscopic  stress-corrosion  crack  is  the  planar  nature 
of  the  crack  surfaces  of  the  glass  fiber  and  matrix  with  little  fiber  pull¬ 
out.  The  time  of  the  brittle  fracture  (final)  stage  of  glass  fiber  and 
surrounding  matrix  is  much  shorter  than  time  tp.  It  follows,  therefore, 
that  the  macroscopic  crack  propagation  rate,  da/ dt,  is  approximately 
given  by: 
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da_D 

dt  tp 


(12) 


where  D  is  the  distance  between  neighbouring  fibers  in  a  doubly 
periodic  square  array  of  fiber  (see  Fig.  4),  written  as: 


D= 


1 


K 


V. 


Substituting  eqns  (11)  and  (13)  into  eqn  (12),  we  obtain: 


(13) 


da 


dt 


1 


sin^  0  f  apF{9)^2Kr  ^ 


V 


Q 


exp 


RT 


Kj 


dO 


J 


) 


(14) 


At  this  point  are  introduced  the  following  quantities: 

4n{  yk  ^ 


V 

J 


ap^27tr 


RT 


so  that  eqn  (14)  can  be  rewritten  as: 


(15) 


da 

dt 


V 


e. 


sin^  0 
0 


ex]p{-iJ,K,F(e)de 


J 


(16) 


Consider  the  result  of  integrating  the  denominator  in  the  above 
equation.  Choosing  sensible  values  oi  (X  =  0.11  ~  0.216  /  mol  (2i), 

Vf=  0.40~0.57,  V  =  0.5  for  a  cross-ply  laminate  and  V  =  1.0  for  a 
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unidirectional  laminate,  R  =  8.31  J/(mol.K)  and  T  -  298  K,  the 
integration  gives  a  value  of  jJ,  estimated  within  the  range: 

jU  =102~521(MPml/2)-i  (17) 

Generally  speaking,  the  range  of  experimentally  observed  stress 
intensity  factor  Ki  at  the  crack  tip  is  between  2  and  15  MPa.mi/2.  it 
follows  that  the  integrand  of  the  denominator  of  eqn  (16)  tends  to  zero, 
except  for  very  small  values  of  0.  Thus,  the  denominator,  which  is 
denoted  by  I,  can  be  written  approximately  as: 

/  =  jg'sin  0 exp (-1.12 Vl  -  cos  d)d6  (18) 


Performing  the  integration: 

r  1  Y  1 


/  =  -4 


A 


l.SSliiK, 


+  0 


exp(-1.58jU0i^;) 


B=e, 


(19) 


-i9=e„ 


0  0 
where  Oq  =  sill-Y?  ~ 


(20) 


Since  the  angle  is  small  and  much  smaller  than  Op,  it  follows  that: 

exp(-1.58^0i,/ir,)  »  exp(-1.58^0,if,)  (21) 


Eqn  (19),  therefore,  reduces  to: 

r  A  \r  I 


/  = 


LV 


l.5^ldKj{h5SldK,  2 


exp(-O.79jU0o^/) 


y 


(22) 


Substituting  eqn  (22)  into  the  denominator  of  eqn  (16),  the  macroscopic 
crack  propagation  rate,  da/dt,  can  be  written  in  the  form: 


dt 


f 


1 


V 


2  +  1.58^^00^/; 


exp(O.79jU0o^/)  (23) 
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This  equation  indicates  that  the  macroscopic  crack  propagation  rate  is 
independent  of  fibre  radius  r  because  /i  is  independent  of  r. 
Furthermore,  the  macroscopic  crack  propagation  rate  is  independent  of 
the  fracture  toughness  of  glass  fiber. 


6.  Results  and  Discussion 
6.1  Verification  by  experiment 

There  is  experimental  data  in  the  literature(i6'  on  a  cross- 
ply  laminate  of  E-glass  fiber  (52.4%  by  vol.)  in  a  diglycidyl  ether  of 
bisphenol-A-type  epoxy  resin  matrix.  The  laminate  consists  of  44  plies 
with  an  individual  ply  thickness  of  about  0.14mm  and  a  stacking 
sequence  of  (0°/90°)iis.  Compact  tension  specimens  with  side  grooves 
were  machined  from  the  laminate.  The  apparent  fracture  toughness  Kc 
of  the  laminate  was  18.2  MPa. mi measured  in  air  at  room 
temperature. 

Stress-corrosion  crack  propagation  tests  were  carried  out  in  0.5N  HCl 
at  298K  under  constant  load.  During  tests,  crack  "mouth"  displacement 
was  measured  at  intervals  using  a  linear  variable  displacement 
transducer.  The  macroscopic  crack  length  was  estimated  from  crack 
mouth  displacement  by  means  of  a  finite  element  method.  Estimates  of 
macroscopic  crack  propagation  rate  were  obtained  by  averaging  the 
derivative  of  the  macroscopic  crack  propagation  length  with  respect  to 
time  over  an  interval  of  between  4  and  24  hours. 

Figure  5  shows  a  logarithmic  plot  of  the  experimental  values  of 
macroscopic  crack  propagation  rate  versus  the  crack  tip  stress  intensity 
factor.  The  figure  reveals  a  linear  log  (da/ dt)  -  log  Ki  relationship  over  a 
an  order  of  magnitude  range  of  da/dt.  Comparison  between 

experiment  and  theory  can  be  made  by  setting  values  of  ^  and  /i,  and 
Oq  and  calculating  the  macroscopic  crack  propagation  rate,  da/ dt,  for 
the  particular  stress  intensity  factor  Kl  using  eqn  (23).  Following 
reference(i6)  where  a  best  fit  method  was  used,  we  take  ^ ,  and  /i  and 
Oq  as  ^=5.01  X  10-15  m.s-1,  jJ,  =  228  (MPa.mi/2)-l  and  Oq  =  0.215°  .  The 
prediction  of  macroscopic  crack  propagation  rate  is  shown  by  the  solid 
line  in  Fig.  5,  and  is  in  good  agreement  with  experimental 
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measurement.  This  gives  confidence  in  the  physical  interpretation  of 
macroscopic  crack  propagation  in  aligned  continuous  GFRP  laminates. 

6.2  Macroscopic  crack  propagation  rate  for  the  limiting  values  of  9q 

The  size  of  a  pre-existing  surface  flaw  on  glass  fiber  affects 
significantly  the  macroscopic  crack  propagation  rate  of  the  laminate. 
Eqn  (23)  indicates  that  the  macroscopic  crack  propagation  rate  increases 

monotonically  with  the  angle  Oq  of  the  pre-existing  surface  flaw  of  the 
glass  fiber. 

First,  consider  the  case  where  Oq  =  0,  where  the  glass  fibers  are  free  of 
surface  flaws.  In  this  case,  eqn  (23)  reduces  simply  to: 

—  =  0.625^p^K^  (24) 

dt 

Substitution  of  eqns  (7)  and  (15)  into  eqn  (24)  yields: 

L25^  ,  (25) 

dt  ' 

Equation  (25)  indicates  the  minimum  macroscopic  crack  propagation 
rate  of  an  aligned  continuous  GFRP  laminate  with  respect  to  Oq,  and 
demonstrates  that  the  variation  of  da/dt  with  Ki  obeys  a  simple  power 
law  having  a  power  of  2.  This  value  is  independent  of  laminate 
geometry  and  acidic  environment. 

Table  1  shows  the  values  of  power  given  in  reported  studies(7^  9,  ii,  12, 15). 
All  the  values  of  power  are  larger  than  2.  It  appears  from  these 
experimental  results  that  the  lowest  value  of  power  is  likely  to  be  2. 

Meanwhile,  let  us  consider  the  case  of  large  values  of  6^  and  consider 
Bartenev's(20)  view  that  the  depth  of  a  pre-existing  surface  flaw  of  fine 
commercial  glass  fiber  is  no  greater  than  0.01  microns.  In  this  case,  the 
angle  0^  would  be  roughly  2.5°  (or  4.36x10*2  rad).  When  6^  is  large 
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and  comparable  with  2.5°,  JlIOqKj  exceeds  unity.  Eqn  (23)  then  reduces 
to: 


^  0-79|/£  exp(O.79Ai0o^/ )  P6) 

In  this  case,  the  macroscopic  crack  propagation  rate,  da/ dt,  is  faster 
than  would  be  expected  using  the  simple  power  law.  The  macroscopic 
crack  propagation  rate  at  high  values  of  stress  intensity  factor  is 
discussed  below. 

6.3  Effect  of  fibre-matrix  interfacial  dehonding  on  macroscopic  crack 

propagation  rate 

For  a  slowly  propagating  macroscopic  crack  there  is  little  or 
no  fiber  pull-out.  In  spite  of  possible  environmental  degradation  of  the 
fiber-matrix  interface,  the  fracture  planes  of  the  fiber  and  matrix  are 
coincidental.  In  contrast,  with  increasing  crack  propagation  rate,  the 
crack  path  is  more  tortuous  due  to  an  increasing  amoimt  of  fiber-matrix 
interfacial  debonding.  The  greater  the  extent  of  fiber-matrix  debonding, 
the  greater  the  likelihood  of  exposure  of  a  larger  surface  flaw  in  glass  to 
environmental  attack. 

Figure  6  shows  the  relationship  between  the  macroscopic  crack 
propagation  rate  and  the  stress  intensity  factor  for  Oq  =  0.2°,  0.3°  and 

0.4°.  The  values  of  ^  and  jd  are  held  constant  as  before.  A  feature  of 
crack  propagation  coexisting  with  fiber-matrix  interfacial  debonding  is 
schematically  illustrated  by  the  broken  line  in  the  figure.  The 
macroscopic  crack  propagation  rate  increases  when  the  stress  intensity 
factor  increases  and  fiber-matrix  interfacial  debonding  occurs  indicated 
by  longer  fiber  pull-out. 

6.4  Effect  of  matrix  toughness  on  macroscopic  crack  propagation  rate 

The  effect  of  matrix  toughness  on  stress-corrosion  cracking 
was  studied  by  Price  and  Hull(i2)  using  the  scanning  electron 
microscope.  They  observed  the  formation  of  matrix  bridges  in  the  crack 
wake  when  tougher,  more  ductile  matrices  were  used.  Such 
mechanisms  that  act  in  the  crack  wake  shield  the  tip  of  the  crack  from 
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local  stress.  Fig.  7  shows  bridging  over  a  crack  of  width  1.  Now  the 
localsed  stress  in  the  matrix  at  the  crack  tip  is  of  the  order  of  it's  yield 
stress.  Consequently,  when  matrix  bridging  occurs,  the  reduced  crack 
tip  stress  intensity  factor  is  given  by(23): 


Ik 


(27) 


By  substituting  eqn  (27)  for  Kl  into  eqn  (23),  the  macroscopic  crack 
propagation  rate  in  terms  of  Ki  is  given  by: 
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dt 


x< 


1 

f 


I  ^ 


2k 


) 


2  +  1. 58^^00 


xexp 


V 
f 

0. 79/1^0  [Kj  -4(7y^^ 


1 


I 


-1 


JJ 


2k 

r\\] 


2k 


J)] 


(28) 


The  numerical  results  of  eqn  (28)  are  presented  graphically  for  1  =  0, 1, 2 
and  3  microns,  Gy^  =  80MPa,  and  the  values  of  /i,  and  6^  were  the 

same  as  those  values  used  previously  (Fig.  8).  The  macroscopic  crack 
propagation  rate  decreases  as  the  size  of  the  matrix  bridge  increases. 
This  bridging  effect  is  more  pronounced  at  lower  values  of  Ki  and  is 
schematically  shown  by  the  broken  line  in  Fig.  8. 


7.  Conclusions 

A  physically-based  micromechanical  theory  of  crack 
propagation  due  to  stress-corrosion  cracking  in  aligned  continuous 
GFRP  laminates  has  been  established.  Good  agreement  between 
experiment  and  theory  give  confidence  in  our  model.  When  the 
relationship  is  expressed  for  the  limiting  case  where  the  glass  fibers  are 
free  of  pre-existing  surface  flaws,  the  relationship  is  represented  by  the 
simple  power  law  with  the  value  of  power  of  2.  The  effects  of  fiber- 
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matrix  interfacial  debonding  and  matrix  toughness  on  the  macroscopic 
crack  propagation  rate  are  taken  into  account  by  the  theory. 


8.  Implications  and  Final  Remarks 

In  the  example  of  predicting  the  behavior  of  a  laminate,  it 
might  be  argued  that  it  would  be  equally  cost-effective  to  simply 
perform  the  mechanical  tests  under  conditions  that  span  the  predictive 
range  of  the  model.  However,  if  the  objective  is  to  design  for  longevity 
and  durability,  then  the  balance  is  shifted  in  favour  of  physical 
modelling.  It  is  in  this  area  that  the  application  of  the  model  is  most 
powerful. 

The  economic  advantage  of  reducing  the  high  cost  of  vast  experimental 
programs  in  assorted  environments  and  stress-states  having  durations 
of  many  thousands  of  hours  is  potentially  huge.  In  addition,  physical 
modelling  provides  the  means  to  assess  the  relative  severity  of  different 
loading  regimes,  e.g.,  constant  amplitude  vs.  variable  amplitude 
spectra,  frequency  effects  and  R-ratio  effects  as  well  as 
load /environment  interactions,  as  described  in  this  report.  This 
capability  is  important  in  the  design  of  experimental  test  programs  that 
ensure  critical  loading  regimes  are  examined. 

Furthermore,  even  though  modelling  of  this  kind  may  be  of  lunited 
benefit  for  predicting,  for  example,  the  notched  fracture  stress  in  quasi¬ 
static  loading  (the  effect  of  lay-up,  hole  size,  environment,  and  so  forth), 
these  models  are  valuable  inasmuch  as  they  draw  attention  to  the 
precise  features  of  the  failure  process  that  have  wider  significance  and 
implications. 

Certain  capabilities  are  missing.  It  does  have  it's  limitations;  in 
particular,  it  lacks  the  ability  to  link  the  observed  damage  growth  of 
splits,  transverse  ply  cracks,  delaminations,  fiber  fractures  growing  in  a 
coupled  manner  to  models  for  the  growth  of  individual  damage  modes 
growing  in  isolation. 

This  modelling  approach,  as  currently  presented,  does  not  have  the 
capability  of  coping  with  variations  in  damage  mode  interactions  from 
system  to  system.  Thus,  whilst  modelling  matrix-dominated  failure 
processes  works  well  for  carbon  fiber  composites  where  the  in  situ 
strength  of  the  0°  ply  in  the  fatigue-damaged  laminate  is  not 
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significantly  different  from  the  strength  of  undamaged  unidirectional 
material,  by  contrast,  an  alternative  model  is  required  for  the  failure  of 
glass  fiber  composites  that  reflects  an  input  for  strength  degradation 
(this  Report). 

There  is  scope  to  integrate  this  modelling  approach  with  the  large 
experimental  programs  that  are  currently  employed  to  design  fracture 
critical  components  and  structures.  However,  this  step  has  yet  to  be 
made.  This  last  point  requires  elaboration.  Modelling  a  particular 
problem  is  only  a  sub-element  of  the  overall  design  process.  I  believe 
that  the  philosophy  behind  the  physical  modelling  approach  has  wider 
applicability.  In  particular,  the  foimdation  of  physical  modelling  could 
be  applied  to  other  problems.  One  important  area  is  based  on 
composite  patch  repair  of  metallic  structures. 

Modelling  the  mechanisms  of  damage  and  relating  them  to 
experimental  data  and  empirical  correlations  that  are  currently 
employed  in  design  provides  a  first  attempt  at  achieving  a  predictive 
capability  in  design,  particularly  with  regard  to  issues  of  strength, 
durability  and  longevity.  It  has  the  potential  to  provide  a  cost-effective 
alternative  to  the  present  largely  empirical  approach  which  takes  an 
excessively  long  time  and  is  expensive.  Even  the  exercise  of  observing 
the  operative  damage  mechanisms  and  calibrating  the  model  for  quasi¬ 
static  loading  provides  useful  physical  insight  which  is  currently  not 
available  to  the  structural  design  engineer  or  materials  technologist. 

However,  the  development  of  this  alternative  design  methodology 
should  be  an  evolutionary  process  with  short  and  long  term  benefits.  In 
the  short  term  current  mechanistic  models  can  be  used  to  bound 
problems,  identify  critical  material  properties,  and  map  likely  failure 
modes.  This  will  lead  to  improved  test  programs  which  concentrate  on 
critical  issues  and  failure  modes  investigated  at  the  important  size 
scale(s). 

Since  the  models  must  be  based  on  experimentally-observed 
phenomena,  experimental  data  will  help  to  improve  and  fine  tune 
them  and  there  will  thus  be  a  synergistic  effect.  In  the  longer  term, 
improved  modelling  will  allow  designers  to  begin  the  design  process 
with  an  intelligent  choice  of  constituent  materials  and  interfaces  and 
then  to  proceed  with  the  prediction  of  performance  and  failure  at  the 
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various  size  scale,  i.e.,  fiber,  ply,  laminate,  structural  component,  full 
scale  structure. 

Existing  design  methodology  at  the  higher  structural  size  scales  will 
thus  be  supported  and  justified  by  fimdamental  imderstanding  at  lower 
size  scales.  Added  benefits  of  this  approach  include  more  options  being 
available  to  the  designer,  a  reduced  need  for  extensive  and  costly 
testing  and  more  efficient  and  shorter  design  iteration  cycles.  This  will 
lead  to  more  versatile,  more  cost-effective,  and  more  efficient 
composite  products,  and  thus  help  to  fulfil  the  promise  of  composite 
materials. 

Finally,  the  key  steps  that  underpin  the  modelling  described  here  are: 
mechanical  testing,  characterization  of  mechanisms  by  observation  and 
modelling  of  these  physical  mechaiusms  to  predict  structural  response. 
If  these  steps  were  routinely  followed  as  part  of  the  design  process  it 
would  potentially  have  a  profound  effect  on  the  overall  cost- 
effectiveness  of  that  process. 
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Fig.  1  Model  of  crack  propagation  in  glass  fibers 
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Fig.  2  Scanning  electron  micrograph  of  fracture  surfaces  of  glass  fibers 
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Fig.  3  Shape  of  stress-corrosion  crack  in  glass  fibers 
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Fig.  4  Glass  fibers  distributed  in  a  doubly  periodic  square  array. 
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Fig.  5  Macroscopic  crack  propagation  rate  versus  the  stress  intensity 

factor  at  the  macroscopic  crack  tip. 
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Fig.  6  Macroscopic  crack  propagation  rate  in  a  state  of  fiber-matrix 
_'i  interface  dehonding. 
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Fig.  7  Matrix  bridging  at  the  crack  front  of  a  macroscopic  stress- 

corrosion  crack. 
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Table  1  Values  of  power  in  a  simple  power  law  equation. 


Lay-up 

Environment 

Value  of  Power 

Ref 

Unidirectional 
E-glass /polyester 

I.ONH2SO4 

3.1 

5 

Unidirectional 
E-glass /polyester 

0.6NHC1 

3.8 

7 

Cross-ply  short 
E-glass /epoxy 

l.ON  H2SO4 

6.3 

9 

Unidirectional 
E-glass  /  polyester 

0.6N  HCl 

3.6 

10 

Cross-ply  short 
E-glass /epoxy 

0.5N  HCl 

3.6 

15 
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Appendix:  Thermal  and  Residual  Stress  Effects  on  Damage 

Accumulation  in  Carbon  Fiber-Epoxi/^  Modelling  the  Failure  Processes 


Damage  Observed  Around  Notches 

An  X-ray  image  of  a  typical  (90/0)s  notched  tensile  specimen  of 
carbon  fibre  in  epoxy  loaded  in  tension  indicates  three  forms  of  matrix 
crack  emanating  from  within  the  notch  tip  damage  zone(i'2): 

(1)  splitting  in  the  0°  plies; 

(2)  delamination  zones  at  90/0  interfaces  whose  size  is  related  to  split 
length;  and 

(3)  transverse  ply  cracking  in  the  90“  phes. 

A  concentration  of  fibre  breaks  in  the  0“  pHes  occurs  at  the  intersection  of 
a  split  and  a  transverse  ply  crack. 

In  a  monotonic  tensile  test  or  in  a  cyclic  loading  experiment,  the  shape 
of  the  damage  zone  remains  unaffected;  it  simply  grows  in  size  in  a 
self-similar  manner,  with  a  characteristic  triangular-shaped 
delamination  zone  at  the  (90/0)  interfaces.  At  the  tip  of  the  split,  the 
delamination  angle,  (between  about  1*  and  7  )  depends  on  ply  thickness, 
on  matrix  toughness,  temperature,  moisture  content.  Essentially,  the 
damage  pattern  at  the  notch  tip  can  be  characterised  by  the 
delamination  angle  and  split  length.  This  ignores  the  transverse  ply 
cracks  in  the  matrix  which  are  important  but  they  can  play  a  secondary 
role. 


Modelling  Damage  in  Monotonic  Loading 

The  change  in  potential  energy  of  a  system  as  a  crack  grows  is 
given  by  the  change  in  strain  energy  less  the  work  done  by  the  applied 
loads.  For  an  infirutesimal  amount  of  crack  growth,  the  energy  released 
xmder  fixed  load  or  fixed  displacement  conditions  is  equal  to: 

5Er  =  -P5u  (1) 

2 

where  P  is  the  remote  load,  u  is  the  imposed  remote  displacement. 
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For  an  increment  of  split  growth,  d  1,  the  energy  absorbed  in  forming 
new  crack  surfaces,  dEab/  is  given  by: 

SEab=Gst5l+Gd(ltana)5l,  (2) 

where  Gs  is  the  energy  absorbed  per  unit  area  of  split;  Gd  is  the  energy 
absorbed  per  unit  area  of  delamination;  t  is  the  thickness  of  a  0  ply, 
and  a  is  the  delamination  angle  at  the  split  tip. 

Some  of  the  global  energy  of  the  system  is  dissipated  when  the  split 
extends  with  a  corresponding  increase  in  specimen  compliance,  dC, 

5Er=^P^5C.  (3) 


For  the  damage  zone  to  grow,  SEj.  >  5Eab_ 

lp^5C  >  Gst5l+  Gd^Sltan  a 

where  P  is  the  applied  load  on  one  quadrant  only 
by 


(4) 


of  the  specimen  given 


P  =  a^(W/2)(2t).  (5) 

Goo  is  the  remote  applied  tensile  stress  on  the  rectangular  tensile 
specimen,  and  W  is  the  specimen  width. 


For  the  limiting  case,  Er  =  Eab, 


(essentially,  G  =  Gc). 


1  P^  5C  ^  ^  Itana 


(6) 


Gc  is  the  effective  toughness  of  the  laminate  and  G  is  the  strain  energy 
release  rate.  A  2-D  finite  element  model  determines  dC/d  i  numerically. 
For  a  range  of  split  length  between  i/a  =  0andi/a  =  6,  it  turns  out 
(fortunately)  that  dC/d i  is  constant,  its  value  dependent  on  the  angle  a, 
the  lay-up  construction  of  the  laminate  and  temperature. 
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Re-arranging  Eqn.  (6),  the  load  for  split  initiation  is  given  by 


2Gst 


il/2 


L8C/8li=oJ 


(7) 


Quite  simply,  the  value  of  Gs  can  be  determined  from  a  measurement  of 
Pi.  For  subsequent  split  growth  under  monotonic  loading; 


p2(5C/5l) 

t  ■ 

ZG^tana 

GdLtana. 

(8) 


Dividing  both  sides  of  Eqn.  (8)  by  notch  length,  a,  and  substituting  for  P 
(in  terms  of  applied  stress  a„),  and  dC/di,  [scaled  from  the  finite 
element  result  {dC/di)\^  ],  this  expression  can  be  normalised  as  follows; 


1  (q..f(t)(W/2)pB(2t)fE  3C 

a  Gdtana(2a/W) 

where  9C  /  3i  is  the  change  in  compliance  of  the  FE  model  with  split 
length;  (W)fe  is  the  width  of  FE  model  specimen;  and  (t)FE  is  the  ply 
thickness  in  the  FE  model. 


Provided  Gs<  Gd  (it  turns  out  that  it  is),  then  the  second  term  on  the 
right-hand  side  of  Eqn.  (9)  is  small,  and  can  be  ignored; 


1  (t)(W/2)pp(2t)i:E  3C| 

a  Gjtana(2a/W)  31 


FE 


(10) 


Eqn.  (10)  indicates  that  for  a  given  2a/W,  the  ratio  1/  a  is  dependent 
only  on  applied  stress,  The  normalised  split  length  with  applied 
stress  curves  for  all  (90/0)s  specimens  should  be  coincidental  regardless 
of  width.  Eqn.  (10)  produces  a  good  fit  to  data(i'2)  using  selected  values 
of  a  =  4*,  Gd  =  400J/m2  and  Gs  =  160  J/m^  for  a  carbon  fibre-epoxy 
composite.  (In  similar  experiments  where  epoxy  was  replaced  with  a 
matrix  of  polyetheretherketone  (PEEK),  a  =  1°,  Gd  =  800  J/m^  and 
Gs  =  400  J/m2).  The  model  can  be  extended  to  fatigue  damage 

accumulation(i'2). 
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Incorporation  of  Residual  Stress  into  the  Finite  Element  Model 

If  thermal  residual  stresses  are  included  in  the  model,  then  a 
simple  comphance  expression  is  no  longer  sufficient. 

As  the  level  of  damage  increases,  three  energy  release  contributions 
exist:  (1)  an  increase  in  mesh-compliance  caused  by  a  change  in  the 
mechanical  strain  energy;  (2)  a  relaxation  of  the  residual  stresses  and 
hence  a  decrease  in  the  thermal  strain  energy  of  the  mesh;  and  (3)  a 
small  expansion  of  the  mesh  due  to  relieving  of  the  residual  stresses  and 
therefore  a  work  contribution  from  the  applied  load.  The  expression  for 
the  total  potential  energy  released  now  becomes; 

dE.  =  -P^9C  -  3Eth  +  P^v  (11) 

2 

Eth  is  the  strain  energy  due  to  residual  stresses  only  and  v  is  the 
displacement  of  the  top  edge  of  the  FE  mesh  before  loading,  measured 
from  the  position  at  i/a  =  0.  At  a  fixed  temperature,  the  elastic 
properties  are  constant  and  Eth  is  proportional  to  the  square  of  the 
remote  residual  stress.  Similarly,  displacement  is  caused  by  damage 
alleviating  thermal  contraction;  therefore  v  is  directly  proportional  to 
the  temperature  drop  applied  and  to  the  magnitude  of  the  residual 
stress.  Incorporating  the  mechanisms  of  splitting  and  delamination: 


ip2  5£_^  +  P^  =  Gst  +  G(j(tana)l  (12) 

2  dl  9l 

The  compliance  change  ac/3i  of  the  FE  mesh  is  independent  of  split 
length.  The  functions  of  thermal  strain  energy  dE^  /  31  and  displacement 
3v  /  3i,  however,  are  not  independent  of  split  length  and  can  be  described 
using  linear  functions  of  1/  a .  The  rate  of  change  of  displacement  is 
independent  of  specimen  size  for  a  constant  2a /W,  but  thermal  strain 
energy  has  to  be  scaled  from  the  FE  result.  Thus: 


r  wt 

3l  31  FE  .WpEtpE 


^  f  lY 

Wt 

WpEtpE. 

(13) 


3v  _  3v 
3l  31fe 


(14) 
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where  C1-C4  are  constants  determined  directly  from  finite  element 
analysis  for  each  temperature  and  material  type.  Substituting  into  Eqn. 
(12)  gives: 


ip2^ 
2  31 


r  ^  r  1 Y 

r  wt 

[WpEtFE. 

+  P 


C3  +  C4 


Gst  +  Gd(tana)l 


(15) 


and  hence: 


^P^(3C/3l)  -  Ggt  -  CiWt  /  Wfe^fe  +  Pf-3 
(G(1  tana  +  C2Wt  /  aW pEtpE  -  PC4  /  a) 


(16) 


If  residual  stresses  are  excluded  from  the  analysis,  C1-C4  are  equal  to 
zero  and  Eqn.  (16)  simplifies  to: 


1=  TP^W3l)-Gst 
G(itana 

Substituting  for  P  and  3C/3i,  and  dividing  both  sides  by  the  notch  half- 
length,  a,  produces: 


1  Yg^Wt(3C/3l)p£WFEtFE  -Ggt-CiWt/  WpEtpE  +<yooWtC3 

a  (Gd  tana  +  C2Wt  /  WpgtpE  —  (JooWtC4) 


(18) 


Calculating  Gs  and  Gd 

Values  of  Gg  can  be  determined  from  the  stress  to  initiate 
damage,  measured  from  x-ray  radiography.  At  the  onset  of 

damage,  the  spht  length  equals  zero  and  Eqn.  (18)  can  be  rearranged  to 
give: 


tTiiW(3C  /  3l)pE^FEtFE 


CiW 

WpE^FE 


+  (T«iWC3 


(19) 


where  substitution  of  appropriate  values  from  finite  element  analysis 
gives  the  energy  required  for  unit  area  of  0“  ply  splitting.  The  work  of 
fracture  for  delamination,  Gd,  is  calculated  for  a  (90/0)s  laminate  at 
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each  temperature  by  fitting  Eqn.  (18)  to  the  experimental  data  for  spht 
growth.  In  some  material  systems  the  delamination  angle  increases 
with  stress  and  a  linear  function  of  the  form  No  has  to  be  substituted  for 
a,  where  N  was  a  constant  determined  for  each  temperature. 

As  the  temperature  is  raised,  the  values  of  Gs  and  Gd  generally 
increase.  The  greater  ductility  of  the  matrix  increases  its  work  of 
fracture  and  the  energy  required  per  unit  area  of  new  crack  surface  rises 
accordingly.  The  values  of  Gd  differ  in  accordance  with  the  toughness 
of  the  matrix  materials. 

The  residual  stress  contribution  is  proportional  to  the  delamination 
angle  a;  a  thin  delamination  zone  gives  rise  to  stiff  coupling  between 
the  0°  and  90°  plies. 


A  Damage  Model  for  Quasi-static  and  Post-fatigue  Strength 

Using  a  finite  element  representation  of  the  damage  at  a  notch 
tip,  combined  with  the  observed  terminal  damage-state  and  measured 
failure  stress,  the  maximum  tensile  stress  in  the  principal  load  bearing 
(0°)  ply  is  approximately  constant  at  laminate  failure,  regardless  of  the 
extent  of  damageO/2),  This  leads  to  a  straightforward  tensile  stress 
failure  criterion; 


rr 


(20) 


where  is  the  remote  failure  stress  of  the  notched  laminate;  Oof  is 
the  localised  tensile  strength  of  the  (0°)  ply  in  the  damage  zone;  and  Kt 
is  the  terminal  notch  tip  stress  concentration  factor.  Oof  is  a  material 
property  which  can  be  determined  independently  while  Kt  can  be 
obtained  from  finite  element  modelling.  For  a  centre-notched  (90/0)s 
specimen  Ktis; 


[l-(2a/wf](Vaf“ 

Furthermore,  the  local  strength  of  a  0°  ply  within  the  notch  tip  damage 
zone  of  a  centre-notched  tensile  (CNT)  specimen  depends  on  the  size  of 
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the  damaged  zone  (just  like  the  strength  of  brittle  solids  depend  on  their 
volume); 


OOf  =  Uq 


Vq 

KVcnt. 


(22) 


where  Gq/  Vq  is  a  mean  failure  stress  and  reference  volume.  KVcnt  Is 
the  equivalent  volume  of  the  0°  ply  in  the  notch  tip  damage  zone,  and  m 
is  the  WeibuU  modulus. 

The  stress  distribution  at  the  notch  tip  varies  with  split  length,  i,  and 
the  delamination  angle,  a,  at  the  tip  of  the  split.  In  my  experiments, 
a  =  3.5%  K  =  0.0014  +  0.0025  ( i/a);  hence 

,.j, _ 

10.0014  +  0.0025(1/  a)8(an)j  1 


=  C6{[0.0014  +  0.0025(l/a)]a^|  (24) 


VcNT  =  8(a2t),  where  t  is  the  thickness  of  the  0°  ply. 

Thermal  strain  mismatch  between  the  0°  and  90°  pHes  causes  a  state  of 
residual  compression  parallel  to  the  fibres  in  the  load-bearing  0  ply. 
This  reduces  the  maximum  tensile  stress  and  hence  reduces  Kt-  As 
temperature  is  increased,  Kt  increases  due  to  two  mechanisms:  (1)  a 
reduction  in  the  level  of  residual  stress;  and  (2)  a  decrease  in  the  lamina 
transverse  and  shear  moduli.  The  latter  causes  an  increase  in  the 
anisotropy  of  the  specimen,  forcing  the  0°  ply  to  carry  a  greater 
proportion  of  the  applied  load.  If  log(Kt)  is  plotted  against  log  (i/a), 
straight  lines  can  be  fitted  to  the  data  for  each  temperature  using  the 
empirical  function  of  the  form: 


^  Cs  (25) 

(Va) 


where  C5  and  p  are  constants. 
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Application  to  the  CNT  Specimen  Notch  Tip 

The  stress  field  at  the  tip  of  a  notch  in  a  CNT  specimen  can  be 
determined  by  finite  element  analysis.  To  determine  KVcnt,  the  finite 
element  results  has  to  be  scaled  from  a  2-layer  quadrant  to  the 
dimensions  of  the  actual  specimen  by  defining  a  reference  volume 
VCNT  =  Sa^to  (  a  is  the  notch  half-length  and  to  is  half  the  total 
thickness  of  the  0”  plies).  K  is  obtained  by  integrating  numerically  using 
finite  element  results  such  that: 


(26) 

i=l 

where  d  is  the  stress  at  a  single  element  integration  point,  Ai  is  the 
corresponding  area  associated  with  the  integration  point,  n  is  the  total 
number  of  integration  points  in  the  finite  element  mesh  and  aOp  is  the 
maximum  stress  in  the  distribution.  KVcnt  therefore  describes  the 
equivalent  volume  of  0“  ply  in  the  CNT  specimen  that  is  subjected  to  the 
maximum  stress  aop. 

As  split  length  increases,  the  stress  gradients  at  the  notch  tip  become  less 
severe  and  a  larger  volume  of  material  is  subjected  to  a  stress  close  to 
the  peak  stress.  An  increase  in  temperature  results  in  a  larger  value  of 
K  due  to  the  increasing  size  of  the  delamination  zone.  A  larger  angle  a 
corresponds  to  more  compliant  shear  coupling  across  the  split  and 
hence  a  larger  equivalent  volume. 

In  order  to  describe  K  as  a  function  of  i/a,  empirical  equations  of  the 
form: 


K  = 


apE 


L^opl 


K  =  C,1- 


(27) 


are  fitted  through  the  data,  where  Ce  is  a  constant. 

If  Eqn.  (27)  is  inserted  into  Eqn.  (22),  the  full  expression  for  the  strength 
of  the  0°  ply  at  each  temperature  becomes: 


-  ^0 


l/m 


C6(Va)8aho 


(28) 
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where  Co,  Vo  and  m  are  measured  from  Weibull  experiments  and  C(,  is 
calculated  using  finite  element  results. 


Notched  Tensile  Strength 

As  subcritical  damage  propagates  from  the  notch  tip,  the  stress 
concentration  factor,  Kt,  continuously  decreases  (Eqn.  (25)).  Damage 
provides  a  progressive  crack  bltmtmg  mechanism  and  the  peak  stress  in 
the  0*  ply  (aop)  is  no  longer  proportional  to  the  applied  stress.  The 
stress  gradients  near  the  notch  tip  are  altered  by  a  change  in  the  stress 
distribution  and  as  the  splits  grow,  a  larger  volume  of  0“  ply  is  subjected 
to  a  stress  near  to  the  peak  stress.  This  in  turn,  causes  a  reduction  of  the 
strength  of  the  0°  ply  (aof)  due  to  a  Weibull  dependence  on  the  volume 
under  stress  (Eqn.  (28)).  Two  competing  mechanisms  therefore  exist  as 
the  splits  extend  -  notch  blvmtmg  and  a  reduction  in  0°  ply  strength. 


Application  of  the  Tensile  Failure  Criterion 

The  remote  stress  on  the  notched  specimen  at  failure  is  given  by: 


V, 

1/m 

■(Va)'' 

.C6(Va)8a\_ 

c,  J 

(29) 


The  remote  failure  stress  is  a  function  of  the  damage  state, 
characterised  by  i/a,  at  the  point  of  failure.  For  typical  values  of  p  and 
m  calculated  previously,  longer  split  lengths  correspond  to  higher 
strength  specimens.  Equation  (29)  therefore  requires  an  experimental 
measurement  of  the  split  length  immediately  prior  to  catastrophic 
failure  to  determine 

For  a  [90/0]s  lay-up,  substituting  for  Va  in  Eqn  (29)  produces: 


tJoof  -  C-ja. 


-2/m 


Wt(9C/9l)pE  WpEtpE  ~  Qst  -  CiWt  /  WpE^FE  + 

(Gdatana  +  C2Wt/ WpEtpE  -aoofWtC4) 


(p-l/m) 


(30) 
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where  C7  is  a  dimensional  constant  which  must  be  determined  for  each 
temperature  and  material  type.  Eqn.  (30)  incorporates  an  implicit 
calculation  of  the  damage  state  at  failure  and  results  in  a  remote 
strength  prediction  which  is  based  entirely  on  independently  determined 
parameters.  The  fracture  energies,  Gg  and  Gd,  are  calculated  by  curve 
fitting  using  the  experimental  data  for  damage  growth,  but  no 
reference  has  been  made  to  the  data  for  tensile  strength. 


Conclusions 

Once  damage  initiates,  the  stress  distribution  at  the  notch  tip  is 
altered  by  the  growth  of  splits  and  delaminations  and  ultimate  fadure  is 
controlled  by  the  maximum  notch  tip  stress  in  the  load-bearing  0  plies. 
Cop.  A  damage-based  model  for  notched  tensile  strength  has  been 
assembled,  using  the  results  of  both  finite  element  analysis  and  Weibull 
experiments.  Thermal  residual  stresses  are  included  into  the  model. 
Additional  energy  contributions  are  suggested  that  may  arise  due  to 
interlaminar  effects  not  included  in  the  2-dimensional  finite  element 
model. 
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